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Abstract 

Common  solutions  for  the  problem  of  random  transient  detection  are  based  on  the 
classical  optimum  detector  for  random  signals  in  noise.  However,  such  a  detector 
does  not  explicitly  use  the  non-stationary  character  of  the  signal  as  a  priori  available 
information.  Reformulation  of  the  optimum  detection  in  the  time-frequency  plane  allows 
one  to  exploit  this  distinguishing  signal  feature  and  suppress  noise.  This  is  accomplished 
here  by  use  of  the  Wigner-Ville  signal  representation  and  an  optimum  signal/noise 
subspace  decomposition  that  improves  the  transient  signal  to  noise  ratio.  The  new 
detection/classification  procedure  eliminates  the  subspace  where  major  part  of  the  energy 
of  random  noise  sample  will  fall. 


1.  INTRODUCTION 

The  problem  of  random  transient  detection  arises  in  the  areas  of  underwater  surveilance, 
seismic  signal  processing,  biomedical  signal  processing,  etc.  Transients  can  be  de¬ 
scribed  as  signals  with  duration  that  is  short  compared  to  the  observation  interval.  They 
can  be  either  deterministic  finite  energy  signals  of  unknown  form,  or  random  finite  dura¬ 
tion  non-stationary  signals.  The  problems  of  detection  and  classification  of  finite  duration 
random  signals,  such  as  the  underwater  acoustic  transients,  have  standard  solutions 
based  on  the  classical  optimum  detectors  of  random  signals  in  noisefl ,  2,  3].  These  de¬ 
tectors  are  not  designed  to  particularly  distinguish  non-stationary  signals  from  stationary 
noise  background.  Here,  we  shall  show  that,  for  the  case  of  Gaussian  transient  detection 
in  Gaussian  noise,  it  is  possible  to  take  advantage  of  this  distinguishing  signal  feature 
by  reformulating  the  optimum  detector  in  the  time-frequency  plane.  This  formulation  is 
based  on  the  Wigner-Ville  signal  representation  and  relies  on  the  generalized  singular 
value  decomposition  to  optimally  separate  the  signal  and  noise  subspaces. 

Classical  optimum  detection  of  Gaussian  random  signals  in  Gaussian  noise  is  reviewed 
in  Section  1 .  Then,  in  Section  2,  the  optimum  detection  procedure  is  reformulated  in  the 
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time  frequency  plane.  The  basic  tool  used  here  is  the  Wigner-Ville  signal  representation. 
Finally,  in  Section  3,  the  generalized  singular  value  decomposition  is  used  to  generate 
an  optimum  stationary/non-stationary  subspace  decomposition.  This  is  the  basis  for  the 
new  detector  which  then  eliminates  the  subspace  where  major  part  of  the  energy  of 
random  noise  sample  will  fall. 


2.  CLASSICAL  OPTIMUM  DETECTION 


We  consider  the  problem  of  detecting  a  zero-mean,  non-stationary,  Gaussian  random 
signal  in  the  zero-mean,  stationary,  colored,  Gaussian  noise  background,  i.e.: 


Ho  :  r  =  n 
H i  :  r  —  ii  -j-  6 


0) 


where  r,  n,  s  are  .Y-dimensional  vectors  of  samples  representing  the  received  waveform, 
random  noise,  and  random  signal,  respectively.  The  noise  is  assumed  to  be  composed 
of  two  independent  components:  colored  interference  and  white  noise  background. 
Therefore,  the  noise  has  a  full-rank  X  x  .V  covariance  matrix  that  can  be  expressed  as 


Rn  =  (t2cRc  +  <r2wI 


where  a2,  are  the  variances  of  the  interference  and  white  noise  components,  respec¬ 
tively.  The  signal  has  '  <  N  covariance  matrix  R3  which  is  possibly  rank-defficient.  It 
is  further  assumed  that  jie  signal  and  noise  are  uncorrelated  and,  therefore,  statistically 
independent.  Then,  the  detection  problem  can  be  rephrased  in  terms  of  the  following 
two  multivariate  normal  models: 


Ho:  r  :  .\  (0.  Rn ) 

Hi:  r  :  Ar(0,  Rn  +  Rs) 


In  words,  the  decision  has  to  be  made  whether  the  received  data  vector  is  a  realization 
of  a  multivariate  zero-mean  normal  distributed  random  vector  with  the  covariance  matrix 

Rn  or  Rn  +  Rs. 


This  problem  has  a  well-known  classical  solution  [1 ,  3]  that  will  be  outlined  next.  We 
shall  then  reformulate  this  solution  in  the  time-frequency  domain,  using  the  Wigner-Ville 
representation,  following  the  approach  in  [4]. 


2.1.  Classical  formulation 

Standard  solution  to  the  detection  problem  above,  based  on  the  likelihood  ratio,  compares 
the  following  sufficient  test  statistic^,  3]: 

7  =  rT\R-{  -(/?„  +  R,)-l}r  >  t  (4) 

Ho 

with  a  threshold  t  determined  according  to  a  suitable  optimization  criterion,  such  as  the 
Ney man- Pearson  criterion,  for  example.  It  is  preferable  to  ’diagonalize’  the  computation 


of  this  statistic  by  means  of  the  generalized  eigendecomposition  of  covariance  matrices 
Rn,  Rs.  Namely,  since  Rn  has  full  rank,  there  exist  a  non-singular  X  x  X  matrix  V' 
such  that[5] 

VtR9V  =  A  =  diag(Xi) 

VT  RnV  =  [  =  diag{  [ }  (5) 

V  =  r{2)  ...  v{S) 


where  A,’s  are  the  generalized  eigenvalues  of  Rs  with  respect  to  Rn,  aranged  in  a  non¬ 
increasing  order  Ai  >  A>  >  . . .  >  Aq  >  0.  Q  =  rank(Rs)  <  X,  and  the  columns  of  V  are 
the  corresponding  generalized  eigenvectors  v{k]’s. 

Now,  a  diagonalizir  j  transformation  can  be  applied  to  the  received  data  vector  r  (with 
dimensionality  reduction,  since  Q  <  X) 

J-  =  VTr  =  [,n  X2  J-3  ...  J-q]T  (6) 

and  the  detection  problem  is  transformed  into  a  decision  which  of  the  following  hypothesis 
is  true 

#o:  .v.X(OA) 

Hi:  x  :  .Y(0, 1  +  A)  U 

The  set  of  generalized  eigenvalues,  {A,},  represent  the  distribution  of  signal-to-noise 
ratio  over  the  coordinates  defined  by  the  generalized  eigenvectors,  SXR  =  trace(A). 
The  test  statistic  becomes 

rj  =  xr[l- (1  +  A)~‘]x  =  jt[A(I  +  A)_1]x  =  xrT.v  >  t  (8) 

J  L  J  H0 

The  resulting  detector  performs  the  linear  transformation,  Eq.  (6),  of  the  received  data, 
with  the  result  being  used  in  the  quadratic  form  computation,  Eq.  (8).  More  explicitly, 
in  terms  of  the  components  of  vector  x  and  the  diagonal  elements  of  I\  A,  the  test 
statistic  is 


EAi  )  * 

i  +  A**7  ~  ^  r 


According  to  Eq.  (9),  the  optimum  receiver  computes  the  O  squared  correlations 
between  the  received  data  and  the  generalized  eigenvectors,  followed  by  their  weighted 
summation.  In  this  sum,  the  signal  components  along  the  coordinates  with  high  signal- 
to-noise  ratio  have  approximately  the  same  weight  %  l,  while  the  components  along  the 
coordinates  with  low  signal-to-noise  ratio  have  reduced  weight  <  1.  The  test  statistic  // 


is  a  quadratic  form  in  Gaussian  variates  with  the  conditional  means  and  variances 


E(r,\H0)  =  trace(T)  =  J2Y^Y 

1=1  1  * 

Q 

E(rj  |  H\)  =  rrace[f(I  +  A)]  =  ^  A, 


Var{r,  |  H0)  =  2  <race(r2)  =  2  j 

Q 

Var(i]  \  Hi)  —  2  trace  T2(I  +  A)"  =  2  ^  A2 

i=i 

The  probabilities  of  detection  and  false  alarm,  PD  and  PpAl  are  governed  by  the 
probability  distribution  of  rj  and  the  threshold  value.  The  distribution  of  the  quadratic 
form  r)  does  not  have  a  neat  analytical  form[2]  and  neither  do  the  probabilities  Pd .  Pfa- 
In  such  a  case,  however,  it  is  common  to  use  the  deflection  signal-to-noise  ratio  as  a 
helpfull  gross  indicator  of  the  detector  performance: 

'Q  X2  j' 

[£(/ 1 //,)-£(/ Iff,))2  [5 

Va,V\H„)  Q  v  ' 

“,=  1  + 

In  the  high  SNR  case,  most  of  the  generalized  eigenvalues  are  large,  A,  »  1,  and 

1  ^  O 

=T  n(S\R,)}‘  »  1  (12) 

When  the  SNR  is  very  low,  a  situation  known  as  a  low  energy  coherence  case  or 
threshold  detection,  all  generalized  eigenvalues  are  small,  A,  <  l,  and 


i  Q  n 

D  *  ^  Ea(  =  v  'lcan  <s-v«r) 


In  summary,  the  classical  optimum  receiver  consists  of  a  bank  of  correlators,  followed 
by  a  weighted  summer  of  their  squared  outputs  (Eq.  (9)),  and  the  threshold  comparator. 
Its  performance  is  governed  by  the  actual  probability  distributions  of  the  quadratic  form  i/ 
under  the  two  hypothesis.  A  gross  indicator  of  the  receiver  performance  is  the  deflection 
signal-to-noise  ratio,  Eq.  (11). 


2.2.  Time-frequency  formulation 

The  optimum  detection  procedure  described  in  the  previous  section  has  an  equivalent 
formulation  in  the  time-frequency  domain.  This  is  of  interest  since  the  signal  to  be 
detected  is  non-stationary,  a  problem  for  which  intuition  suggests  comparison  of  time- 
frequency  signal  representations^]. 

The  time-frequency  domain  equivalents  of  the  inner  products  in  the  classical  optimum 
receiver,  Eq.  (9),  can  be  obtained  using  any  bilinear  transformation  that  preserves  the 
inner  products.  Among  such  transformations,  Wigner-Ville  representation^]  has  certain 
advantages  that  make  it  the  most  attractive  for  our  purposes.  Namely,  it  is  the  only  one 
that  has  all  of  the  following  desirable  properties[6,  7]:  (1)  it  is  real-valued;  (2)  it  is  invariant 
to  time  and  frequency  shifts;  (3)  it  has  the  same  region  of  support  as  the  signal;  (4)  it 
preserves  the  inner  products;  (5)  it  has  the  proper  marginals;  (6)  among  all  the  bilinear 
representations  satisfying  (1)— (5)  it  is  the  one  best  localized  in  the  time-frequency  plane. 

It  should  be  noted,  however,  that  these  properties,  while  valid  in  the  continuous-time 
case,  are  not  generally  valid  for  the  discrete-time  Wigner-Ville  representation  which  is 
periodic  in  frequency  with  period  s-  [6].  All  of  the  above  properties  are  valid  only  if  the 
discrete-time  signal  spectrum  is  non-zero  over  an  interval  less  than  -k  on  its  fundamental 
period  2k  [6].  Two  practically  important  such  cases  are:  (a)  when  the  real-valued 
continous-time  signal  is  oversampled  by  at  least  a  factor  of  2,  and  (b)  when  the  signal  is 
analytic,  i.e.  its  spectrum  vanishes  for  the  'negative’  frequency  half  of  the  fundamental 
period  2 7r  [6].  In  the  rest  of  this  paper,  it  will  be  assumed  that  the  signals  of  interest 
satisfy  the  aforementioned  condition. 

For  signals  of  interest  that  have  a  finite  discrete-time  support  region  T  of  length  .V  = 
2/v  +  l,  i.e.x{k)  =  y(k)  =  0.  k  g  T,  the  discrete  Wigner-Ville  representation  is  given  by 

K 

Wx{k,  l)  =  2  J2  x(k  +  m)  J*( k'  ~  m)  (1 4) 

m=  —  h 

where  k  e  T,  and  /  e  fl  represents  M  >  X  frequency  domain  sampling  points. 

In  the  absence  of  aliasing,  the  following  discrete-time  discrete-frequency  version  of  the 
Moyal’s  formula  is  obtained 

-  (jf  ••'•')('?  •'<)*  (l5> 

1  len  *eT 

For  example,  this  form  of  the  Moyal’s  formula  is  applicable  to  real-valued  signals  over¬ 
sampled  by  at  least  a  factor  of  two,  and  to  analytic  signals.  This  is  the  equation  that  allows 
easy  reformulation  of  the  classical  optimum  detection  procedure  in  the  time-frequency 
plane.  However,  it  is  necessary  to  consider  the  Wigner-Ville  representation  of  random 
signals[8]  first. 

Let  x  be  a  realization  of  a  harmonizable  zero- mean  complex-valued  discrete-time  random 
process  X  with  finite  discrete-time  region  of  support  T  of  length  .V  =  2A  +  1,  e.g.  finite 


length  random  transient.  Then,  for  k  e  T 

Rt{k  +  m,  k  —  m)  =0,  \m\  >  K 

and  the  following  expectation  exists[8] 

K 

W.v(fc,  /)  =  E{  Wx(k,  /)}  =  2  E{x(k  +  m)  x*{k  —  m)}  e 

m=—K 

K 

=  2  Rz(k  +  m,  k  —  m)  , 


-J2~hn 

k  €  T 
i  e  ft 


and  is  called  the  Wigner-Ville  spectrum  of  the  process  .  Here,  T.  ft  are  the  discrete¬ 
time  and  discrete-frequency  support  regions,  respectively. 

Infinite  duration  random  processes  are  usually  observed  over  finite  intervals,  and  their 
covariances  may  be  known  or  estimated  only  on  a  grid  of  limited  extent.  In  such  a  case, 
it  is  possible  to  obtain  an  approximation  of  the  sampled  Wigner-Ville  spectrum,  smeared 
one-dimensionally  in  the  direction  of  frequency,  by  means  of  the  discrete  Wigner-Ville 
spectrum,  Eq.  (17). 

For  example,  let  X  be  a  harmonizable  infinite  duration  discrete-time  random  process 
observed  on  the  interval  -IK  <  k  <  2I\,  with  spectral  support  less  then  -  and  the 
covariance  Rz(p.q)  known  on  the  grid  -2  A'  <  p.q  <  2K.  Then,  Rx(k  +  m,k~m)  is 
completely  known  on  the  grid  -I<  <  k.m  <  I<  and  the  sampled  approximation  of  the 
Wigner-Ville  spectrum  can  be  computed  as 


W.vik.l)  =  2  g(m)  Rz(k  +  m.  k  -  in)  e~J%m 


=  2  £  Rx(k  +  m,k-  m)  lefi.  -K  <  k  <  K 

m=  —  K 

where  g(m)  is  the  rectangular  truncating  window  with  Ci(l)  as  its  discrete  Fourier  trans¬ 
form.  If  the  process  X  is  also  stationary,  since  there  is  no  aliasing  in  the  time-frequency 
plane,  the  Wigner-Ville  spectrum  equals  the  ordinary  power  spectrum,  and  its  sampled 
smeared  approximation  is 

Wx(k.l)  =  G(l)  *  Sz(l)  =  Sz(l),  (19) 

In  words,  the  finite-record  approximation  of  the  Wigner-Ville  spectrum  of  an  adequately 
oversampled  stationary  process  is  equal  to  the  finite-record  approximation  of  its  ordinary 
power  spectral  density. 

Now,  the  detection  problem,  Eq.  (1),  can  be  rephrased  in  the  time-frequency  plane  as 
the  choice  between  the  hypotheses 

"»  :  ,,V  -  ll"  (20, 

h,  :  wr  =  u„  +  '  ’ 


where  vr„,  IF,  are  the  discrete  Wigner-Ville  representations  of  the  particular  noise  and 
signal  realizations. 

We  are  ready  now  to  transform  the  time-domain  correlator  based  optimum  receiver,  Eq. 
(9),  into  the  time-frequency  domain  correlator  based  receiver.  The  inner  products  in  Eq. 
(9)  can  be  replaced,  according  to  the  discrete  Moyal’s  formula,  Eq.  (15),  by 

rT  ■  <->■']'  =  ^TjY.Y.  lV '(*■■')  »'.(*■■ ')  <21> 

“*u  ten  ker 

where  Wr(k.  /)  is  the  Wigner-Ville  representation  of  the  received  data  vector  and  \Yt(k.  I ) 
are  the  Wigner-Ville  representations  of  the  generalized  eigenvectors  rl,).  Using  this 
result  in  Eq.  (9),  we  get 

=  kJ) H>< f  (22) 
lea  ke t 

where  B(kJ)  is  the  weighted  sum 


1  +  A, 


W,(k.l) 


which  is  fixed  by  the  generalized  eigendecomposition,  Eq.  (5),  does  not  depend 
on  the  data  vector  r,  and  can  be  precomputed.  If  the  slices  of  \Vr  and  B  for  the 
fixed  time-index  k  are  labeled  wj  =  [U'r  ( A*.  0)  \Yr(kA)  ...  \Vr(k.  M  -  l )]  and  b[  = 
[B(k.  0)  B(k.  i)  ...  B(k.  M  -  I)],  respectively,  the  optimum  test  has  a  more  explicit 
form 

n  =  l>[  ■  >rk  >  /  (24) 

*1t 

Accordingly,  the  time-frequency  domain  equivalent  of  the  classical  optimum  receiver 
computes  the  slice  by  slice  cross-correlation  of  the  received  data  Wigner-Ville  repre¬ 
sentation  with  the  precomputed  weighted  sum  of  Wigner-Ville  representations  of  the 
covariance  matrix  generalized  eigenvectors.  An  alternative,  expression  for  the  optimum 
test  is  obtained  by  introducing  the  matrix  notation  U'r  =  [U'rU\/)].  B  =  [B(k.l) ] 

rj  =  trace  l  B  Wf  \  —  tracrl  Br  Wr\  J  t  (25) 


In  the  reference  [4],  a  similar  reformulation  of  the  continuous-time  optimum  detector  has 
been  presented  for  the  case  of  white  noise  only. 

The  test  statistic  conditional  expectations  remain  the  same  as  those  of  the  classical 
receiver,  as  given  by  Eqs.  (10).  It  is  interesting  to  relate  those  results  with  the  Wigner- 
Ville  spectra  of  signal  and  noise.  Taking  the  conditional  expectations  of  Eq.  (25),  using 


the  Eqs.  (17)  and  (18),  and  combining  the  results  with  Eqs.  (10)  we  get 


E(n  I  Ho)  =  tracers  Wj)  =  £ 

i=l 

_  Q 

E(rj\H i )  =  trace  [tf  (VV^  +  Wjj]  =  A, 


Var(r)  \  Hn)  =  2  trace 


1=1 

Q 


(bw5)%2x:(t^- 


(26) 


\  'ar(r]  |  Hi )  =  2  trace  {[*  (vvj  +  wj)]'|  =2X;Al2 


Q 


where  Hy  is  the  finite-record  approximation  of  the  Wigner-Ville  spectrum  of  noise,  as 
discussed  in  connection  with  Eq.  (18),  and  Ws  is  the  discrete  Wigner-Ville  spectrum  of 
a  finite  duration  non-stationary  signal,  as  in  Eq.  (17).  The  deflection  signal-to-noise  ratio 
is  the  same  as  before  and  given  by  Eqs.  (11)  —  (13). 

We  see  that  an  alternative  time-frequency  domain  realization  of  the  optimum  receiver 
is  possible.  It  matches  optimality  with  physical  interpretation  provided  by  the  use 
of  the  Wigner-Ville  representation,  which  is  particularly  relevant  for  the  detection  of 
non-stationary  signals.  However,  the  receiver  performance  remains  the  same.  The 
real  significance  of  this  time-frequency  formulation  is  that  it  will  allow  incorporation  of 
additional  information  that  may  bring  improvement  in  performance,  as  will  be  shown  in 
the  next  section. 


3.  SUBSPACE  BASED  DETECTION 

Our  objective  is  to  use  a  priori  available  information  that  signal,  as  opposed  to  the  noise, 
is  non-stationary  in  order  to  suppress  that  background  noise  and  thereby  improve  the 
deflection  signal-to-noise  ratio  and  detector  performance.  How  do  we  accomplish  this 
is  discussed  next. 

3.1.  Optimum  signal/noise  subspace  decomposition 

The  Wigner-Ville  spectrum  of  the  stationary  random  signal  and  its  finite  record  approxi¬ 
mation,  Eq.  (19),  are  simple  outer  products  of  the  power  spectral  density  as  a  function 
of  frequency  and  a  constant  function  of  time.  Consequently,  the  matrix  of  samples  of 
the  finite  record  approximation  of  the  Wigner-Ville  spectrum  of  stationary  noise,  wy,  is 
a  unit-rank  matrix.  On  the  other  hand,  the  discrete  Wigner-Ville  spectrum  of  the  random 
transient,  Eq.  (1 7),  represented  by  the  matrix  Wy  has  rank  that  is  generally  much  higher 
than  one.  This  distinction  can  be  exploited  to  suppress  stationary  noise  and  improve  the 
detector  performance. 


The  generalized  singular  value  decomposition  of  the  signal  Wigner-Ville  spectrum  with 
respect  to  the  noise  Wigner-Viile  spectrum  provides  a  necessary  tool  to  accomplish  this. 
Assuming  that  Vty  and  VV5  are  V  x  A/,  A/  <  .V  matrices,  then  there  exist  a  non-singular 
.1/  x  M  matrix  X  and  orthogonal  matrices  Us(\  x  .V).  r„(.V  x  X)  such  that[5,  9] 

VV5  =  USDSX~{ 

vvy  =  cnDnx~' 

Ds{ A'  x  A/)  =  diag(at).  0  <  c>i  <  o_>  <  ...  <  01  /  (27) 

£)„(A'  x  .V/)  =  diag(3, ),  >  ...  >  3q  >  =  ...  =  A.i/  =  0. 

+  J,2  =  1,  <7  =  rank{W\') 


where  j(,),  the  columns  of  A',  are  the  generalized  singular  vectors  of  wy  and  VV5 
satisfying 

3;  wj Ws  -r("  =  a2  VV^W.v  -c(i)  (28) 

Note  that  for  the  stationary  noise  rank- (wy)  =  q  =  1  and,  therefore,  only  y  y  0. 

Let’s  introduce  the  following  matrix  partitions 


.v  =  [j-111  A',].  11  (•„]. 

(•„,]■  d,  =  [n0' 


(29) 


where 


Then, 


A’i  =  Ai(  A/  x  M  -  1). 

D\  =  D\(.\  —  1  x  A/  —  1 )  =  diag{a-2 . o.v) 

=  r,i(x  x  .v  - 1 ). 

rni  =  r„,( a  x.v  -i) 

WsX  =  CsDs  =  [aiu'"  ra\D\] 

W'.V'A  =  l’nDn  =  (),] 


where  Oi  is  a  A  x  M  -  l  zero  matrix.  More  specifically, 


VVsX, 

W.v'A  1  =  O, 


(30) 

(31) 

(32) 


i.e.,  A'i—  transformation  maps  the  Wigner-Ville  spectrum  of  stationary  colored  noise  into 
an  all-zero  matrix  Ot.  At  the  same  time,  since  a]  <  «,  is  small,  the  signal  Wigner-Ville 
spectrum  is  mapped  with  the  minimum  loss  of  signal  power. 

We  can  interpret  these  results  as  saying  that  all  of  the  noise  energy,  together  with  a 
minimal  fraction  of  signal  energy,  lie  along  the  first  generalized  singular  vector 
and  there  is  no  noise  energy  in  the  subspace  determined  by  the  remaining  generalized 
singular  vectors  that  constitute  the  columns  of  A'i[9].  Therefore,  we  can  suppress  noise 
by  projecting  the  Wigner-Ville  spectrum  onto  the  subspace  spanned  by  the  columns  of 


A'i-  To  accomplish  this,  we  need  the  subspace  projection  matrix  that  can  be  computed 
by  the  following  stable  procedure. 

Quite  often,  the  matrix  X-1  is  ill-conditioned  and  its  inversion  produces  inaccurate 
generalized  singular  vectors.  Since  we  are  mainly  interested  here  in  the  subspaces 
spanned  by  the  generalized  singular  vectors,  numerically  more  stable  procedure[5], 
adequate  for  our  purposes,  is  to  use  the  QR  decomposition  of  this  matrix 

X~T  =  QR.  X~l  =  RtQt  (33) 

Here  R  is  an  upper  triangular  matrix,  and  Q  is  an  orthogonal  matrix  whose  first  column 
is  the  unit  vector  in  the  direction  of  the  first  generalized  singular  vector  while  the  remain¬ 
ing  columns  define  an  orthonormal  basis  for  the  subspace  spanned  by  the  remaining 
generalized  singular  vectors  x{i).  i  =  2.3 . M.  The  generalized  singular  value  de¬ 

composition,  Eq.  (31),  now  takes  the  form 


WsQ  =  i'sDs  Rt 

=  l'nDn  RT 

If  we  use  the  partition 

Q  =  W"  Qi) 

where  qll)  is  the  first  column  of  Q  while  Q  t  contains  the  remaining  ones,  then 

P  =  QiQj  (36) 

is  the  projection  matrix  onto  the  subspace  spanned  by  the  columns  of  Xj.  Then,  the 
projection  of  the  noise  Wigner-Ville  spectrum  onto  this  subspace  equals  zero, 

VV\-P  =  0  (37) 

Therefore,  using  this  projection  operator  the  Wigner-Ville  spectrum  of  the  stationary  noise 
is  anihilated  while  that  of  the  transient  signal  is  barely  affected. 

3.2.  Detection  with  noise  subspace  elimination 

The  generalized  singular  value  decomposition  of  WA-  and  Ws,  Eq  (27),  thus  generates 
an  optimum  stationary/non- stationary  subspace  decomposition,  Eq  (31),  tha*  will  be  the 
basis  of  the  new  detection  scheme  discussed  next. 

Rather  than  correlating  the  received  data  Wigner-Ville  representation  with  a  reference,  as 
in  Eq.  (25),  we  propose  to  first  use  the  noise  suppressing  projection  onto  the  subspace 
spanned  by  the  columns  of  Xt .  It  is  easy  to  show  that 

H0:  E{\VrQl)  =  V^vQ,  =  0] 

W,  :  E[\YrQi)  =  VVv'Q,  +  WsQi  =  Y 


(34) 

(35) 


(38) 


Therefore,  we  propose  to  compute  the  test  statistic 

n  =  ^•«ce|(5(?i)(lir^i)r| 

=  trace  j  B\  H'/'l  >  t 1 

t  i  H0 


(39) 


where  tx  is  the  threshold  determined  by  a  suitable  criterion,  while  using  P,  the  projection 
matrix  onto  the  subspace  spanned  by  the  columns  of  AT  (Eq.  (36)),  we  have 

B i  =  B{QiQj)  =  BP  (40) 


The  conditional  means  can  be  expressed  as 


E(i)\Hq)  =  trace 
E{i]\H\)  =  trace 


BiW’Xj  =  0 

B  i  ( VV  v'  +  VV5)rj 


=  trace  (  B\  W'J  j 


(41) 


However,  the  conditional  variances  cannot  be  brought  into  a  neat  form,  but,  intuitively 
we  expect  them  to  be  reduced  with  respect  to  the  classical  case  due  to  the  elimination 
of  noise  subspace  by  the  projection  onto  the  AT— subspace.  This  point  is  verified  by 
Monte-Carlo  simulations  discussed  next. 


3.3.  An  example 

As  an  illustration,  we  consider  a  case  where  signal  and  noise  are  specified  according  to 
Fig.  1 .  The  signal  autocovariance  and  Wigner-Ville  spectrum  are  given  in  Figs.  1  .a  and 
l.b,  respectively.  Power  spectra  of  the  signal,  noise,  and  signal  plus  noise  are  given  in 
Figs.  1  .c,  1  .d,  and  1  .e,  while  the  Wigner-Ville  spectrum  of  the  noise  is  shown  in  Fig.  1  .f. 
Narrowband  interference  is  five  times  stronger  than  the  white  noise  component.  Overall 
signal  to  noise  ratio  is  -7.8dB,  while  signal  to  white  noise  ratio  is  OdB.  Eigenspectra 
of  the  signal  and  noise  covariance  matrices  are  shown  in  Figs.  2.a  and  2.b.  Rank  of 
the  signal  covariance  is  14,  while  the  noise  covariance  has  full  rank  32.  Singular  value 
spectra  of  the  signal  and  noise  Wigner-Ville  spectra  (WVS)  are  presented  in  Figs.  2.c 
and  2.d.  Rank  of  the  signal  WVS  is  31  and  that  of  the  noise  WVS  is  1.  Generalized 
singular  values  of  signal  and  noise  WVS’s  and  generalized  eigenvalues  of  signal  and 
noise  covariances  are  shown  in  Figs.  2.e  and  2.f.  Finally,  the  corresponding  WVS 
matching  function  B,  Eq.  (23),  and  the  WVS  subspace  matching  function  Bu  Eq-  (40), 
are  displayed  in  Figs.  3. a  and  3.b.  Alternatively,  the  effect  of  projecting  (filtering)  the 
signal  and  noise  WVS’s  onto  the  no-noise  subspace  is  shown  in  Figs.  3.c  and  3.d. 

In  order  to  verify  the  conjecture  made  at  the  end  of  the  previous  section,  the  operation  of 
both  the  classical  and  the  new  detector  on  the  sample  vectors  of  signal  and  noise  defined 
above  was  simulated  with  500  Monte-Carlo  runs.  While  this  number  of  simulations  is 
insufficient  for  accurate  estimates  of  the  probabilities  of  detection  and  false  alarm,  it  was 
sufficient  to  estimate  the  deflection  signal  to  noise  ratios  of  two  detectors.  The  estimated 


deflection  signal  to  noise  ratio  gain  obtained  in  this  manner  was  1.9  (or  about  3dB), 
confirming  our  conjecture  that  the  new  detector  provides  additional  noise  suppression. 
Of  course,  this  is  far  from  being  a  definite  proof  of  superiority  of  the  new  detector  under 
all  conditions.  A  more  comprehensive  set  of  simulations  is  necessary  for  more  definite 
statement  of  such  kind.  It  should  be  emphasized,  however,  that  the  noise  suppression 
procedure  described  here  is  non-parametric,  relies  only  on  noise  stationarity  feature, 
and  therefore  should  be  quite  robust  providing  gain  under  various  different  signal/noise 
scenarios. 


4.  SUMMARY 

A  new  procedure  for  detection  of  Gaussian  transients  in  stationary  colored  Gaussian 
noise  was  presented  that  exploits  the  signal  non-stationarity  as  its  distinguishing  feature. 
The  classical  optimum  detector  for  Gaussian  random  signal  was  first  reformulated  in  the 
time-frequency  domain  using  the  Wigner-Ville  representation.  In  this  domain  the  distinc¬ 
tion  between  non-stationary  and  stationary  processes  is  easy  to  exploit.  The  generalized 
singular  value  decomposition  of  the  signal  and  noise  Wigner-Ville  spectra  generated  the 
optimum  signal/noise  subspace  decomposition.  By  an  appropriate  subspace  mapping, 
the  noise  subspace  was  eliminated  providing  an  extra  degree  of  noise  suppression.  This 
point  was  verified  by  simulations  for  a  particular  set  of  conditions.  Although  the  new  pro¬ 
cedure  appears  to  be  robust  with  respect  to  changing  signal/noise  situations,  a  more 
comprehensive  set  of  simulations  is  needed  for  a  definite  conclusion.  The  approach  ad¬ 
vanced  here  can  easily  be  extended  to  the  optimal  linear-quadratic  detector  for  problems 
involving  certain  classes  of  non-Gaussian  background  noise[10]. 
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(b)  WVS  Subspace  Matching  Function 
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